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Abstract. We analyze the crossing of a quantum critical point based on exact 
results for the transverse XY model. In dependence of the change rate of the 
driving field, the evolution of the ground state is studied wrhile the transverse 
magnetic field is tuned through the critical point with a linear ramping. The 
excitation probability is obtained exactly and is compared to previous studies and 
to the Landau-Zener formula, a long time solution for non-adiabatic transitions 
in two-level systems. The exact time dependence of the excitations density in 
the system allows to identify the adiabatic and diabatic regions during the sweep 
and to study the mesoscopic fluctuations of the excitations. The effect of white 
noise is investigated, where the critical point transmutes into a non-hermitian 
"degenerate region" . Besides an overall increase of the excitations during and at 
the end of the sweep, the most destructive effect of the noise is the decay of the 
state purity that is enhanced by the passage through the degenerate region. 
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1. Introduction 

Adiabatic passages are a fascinating and important tool in modern physics, since they 
aUow to manipulate quantum states in a controlled manner via tunable parameters 
of a physical model system. Their potential impact ranges from setting current 
standards by means of adiabatically pumping charge through a quantum dot [H [2l [3l [4] 
to the access to ground states of e.g. non-integrable Hamiltonians via adiabatically 
connecting to it from an integrable Hamiltonian[5]. Since modern methods in quantum 
optics allow fabrication of a wide class of lattice Hamiltonians with systems of 
atoms loaded into a suitably engineered optical trap, these techniques have become 
experimentally accessible nowadays. We here will discuss the latter application 
of a possibly adiabatic transit of a ground state with an explicit time dependent 
Hamiltonian and "slowly" varying parameters, which in the recent literature has 
also been termed adiabatic quantum computation (AQC)[5]. The efficiency of this 
procedure is compromised by size effects and by the presence of noise, the arch-enemy 
of all quantum information tasks. 

The closing of the gap in the spectrum of excitations with growing system size 
gives a limit to the rate the parameters may be changed in order to have a sufficiently 
high fidelity with the ground state of the final Hamiltonian. In other words, the 
gap rules the time of computation. An intriguing question is then, what happened 
if a quantum critical point [6j is crossed during the passage. In the dynamic scenario 
presented here, the vanishing gap makes the notion of slow change rates obsolete, 
and strictly speaking no adiabatic passage exists across a quantum critical point. In 
the laboratory we have to deal with finite systems and as a consequence the gap, 
though scaling down to zero somehow with the system size, will nevertheless be 
finite. Then, the question arises how the fidelity of the ground state and/or the time 
of computation scales with the system size. This question is relevant for deciding 
whether adiabatic quantum computation may be successful or not. It has been 
addressed for the one-dimensional transverse Ising and the anisotropic Heisenberg- 
chain[8] in numerically solving the Heisenberg equation; the excitation probability has 
been calculated perturbatively taking into account a finite number of multiparticle 
excitations. The authors concluded that the adiabatic time scale for the Ising model 
in a transverse field scales polynomially with the system size. Besides a complete 
treatment of many-particle excitations this work did not take into consideration the 
crucial role fiuctuations of the average excitation number play in particular for small 
system size. The previous question is also related to the Kibble-Zurek theory, where 
critical scaling laws have been employed to obtain predictions for the density of defects 
after the crossing of a real phase transition[TDl[TTl[T2]. Its predictions have been tested 
for the one-dimensional transverse Ising model in Ref. [9] , where the kink density has 
been calculated numerically and good qualitative agreement with the Kibble-Zurek 
scenario has been reported up to little quantitative discrepancy. By means of simple 
scaling arguments the density of created excitations as function of the sweeping rate 
has been evaluated for models belonging to different universality classes [TB]. Using 
the exact solution for the transverse Ising model [7l [Ml [151 [IF' j the kink density at 
vanishing magnetic field has been calculated in Ref. [17j , and numerical discrepancies 
with Ref. :9, have been noted, where the time dependent Bogoliubov approach has 
been applied and excitation probabilities in the long time limit are obtained from a 
mapping to Landau-Zener tunneling. 

How the presence of the noise affects this scenario is a rather unexplored, though 
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very important, topic. In experimental setups one has to struggle with the influence of 
environmental degrees of freedom. This transforms the avoided criticality of the finite 
system into a region of non-hermitian degeneracy [H] extending over a finite interval in 
control parameter space. This means that the system passes through a gapless phase 
within a finite time window. Moreover, besides the inverse of the spectral gap that 
sets a lower bound on the sweeping period, the noise introduces a second time scale: 
the dephasing time, that on the contrary sets an upper bound. These competing time 
scales give rise to non-trivial effects, which are largely ignored in the recent literature 
on AQC [ini HOI m] notwithstanding their importance for a realistic analysis of AQC 
protocols. 

A prominent and simple class of models with quantum critical point are the one 
dimensional XY models in transverse magnetic field[6l [7]. In this work, we trace 
the exact evolution of the model during a sweep in order to study the robustness 
of the adiabatic passage against size effects with and without the presence of noise 
in the driving field. The paper is organized as follows: In section [5] we present the 
model Hamiltonian together with a sketch of its exact solution of the dynamics and we 
briefly review the Kibble-Zurek scenario for a quantum phase transition in section [3l 
In section|4j the results for the excitation probability and its mean square fluctuations 
during the sweep are presented and discussed. In section [5l we analyze the case of a 
noisy driving field. Conclusions out of our findings are drawn in section [B] 



2. The model and its exact solution 



The model Hamiltonian under consideration throughout this work is the one- 
dimensional spin- 1/2 XY model in transverse magnetic field, an archetype model 
exhibiting a quantum phase transition[6 for which an exact solution exists [71 1141 

MM 



H 



L 

E 



(l + 7)5J5jVi + (l-7)^j'^jVi 



h{t) y: 



(1) 



^From here on, for sake of simplicity the unit of energy, h and the lattice spacing 
are set to one. In the above equation 7 is the in-plane anisotropy coefficient, 
the magnetic field h{t) is the time-dependent tunable parameter and the boundary 
conditions are periodic. The model undergoes a quantum phase transition at — I 
separating a paramagnetic {h > he) from a broken symmetry phase, where (S^) ^ 0. 
Details about the exact solution and the derivation of the exact time evolution 
are found in [7], and the main steps are outlined in the Appendix. The canonical 
transformations (|A.1[) . (jA.4p and (|A.7p completely decouple the Hamiltonian into a 
direct sum H = (S)\k\Hk of 4-dimensional Hamiltonians Hk acting non-trivially only 
in the Hilbert space {|0)^ _f. , |fc, —k) ; |A:) , |— A;)}, where, due to the parity symmetry, 
it further decouples into Hk = H^'^'^ ® H^^*^^ . Whereas in the odd-occupation Hilbert 
space the Hamiltonian is already diagonal, H^'^'^ = — cos0fel2, in the even-occupation 
Hilbert space and basis {|0) , |fc, — fc)} one obtains 



Tjeven 



■ C0S(?;>fcll2 -t- flfcCT^ + 5fcCr^ 



(2) 

where (^k — i'Kk/^L — 1), ak = cos (j)k + h{t)^ bk = 7sin(^fc, and cr" are the Pauli 
matrices. /.From now on we will focus on the nontrivial even part of the Hamiltonian 
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and for sake of simplicity the superscript "even" will be dropped. The eigenvalues are 



E,^±Ak^±^ai + bi (3) 
with Ok and bk as defined above, and the ground state is 

\GSk) - hf") =: i-ifh + |0),._, , (4) 



where 



Afc + Ofc 

Q^fe = .... , Pfc 



and the fermionic operators // are defined in Ea. (|A.7p . The ground state of the whole 
chain, which later on will be the initial state, is the tensor product 

® <» 
\GS) = H \GSk) , ft.o = n \GSk) {GSk\ , 

\k\ \k\ 

since we will start in the phase where the parity symmetry is not broken. As a 
consequence, (S^) — during the time evolution. 

In terms of single mode operators the Hamiltonian can be written as 

H = J2Ak{vlvk-l) , (5) 

fc 

where the relation between the fermionic operators rjk and fk reads 

Vk = /3k fk + iakflk ■ 
It is worth noticing that the diagonalization is not giving the time evolution, since 
[dtH, H] ^ 0. The time evolution operator is given by 

U{t) ■.= Texp{-i [ dTH{T)} (6) 
Jo 

with T the time ordering operator, and we have to deal with the Schrodinger equation 

® 

idtU^HU] [/ = [|[/('=^ (7) 
|fc| 

for the time evolution operator U instead. An exact solution does exist for several time 
dependences [7l[22j[23]. In the simplest case of Hnear sweep h{t) = ho + {hi — ho)t/T 
the long time physics exhibits the Landau-Zener tunneling scenario [221I21]; for details 
of the exact time evolution we refer to the Appendix. 

The fidelity of an adiabatic passage process is the overlap of the wavefunction 
at time t and the ground state \GS{t)) of H{t) at the same time. The latter would 
be the wavefunction for adiabatically driven magnetic field. The overlap integral 
{GS{t) \ U{t) \GS{0)) is thus a rough measure for the efficiency of the adiabatic passage. 
Another way to check the adiabaticity of a sweeping process is to evaluate the average 
number of excitations in the system after the sweep, i.e. the expectation value of 
J2k VkVk on the state U{t) \GS{0)). The passage is fully adiabatic if no excitations are 
produced at the end of the drive. 
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3. Kibble-Zurek scenario 



The density of excitations crossing a phase transition can be foreseen according to 
the Kibble-Zurek scenario [lOl [TTl [12] . Before entering the bulk of our work, we shall 
give a short summary of results from the recent paper by Zurek, Dorner and ZoUer 
(ZDZ) [9j, who extended the Kibble-Zurek scenario to quantum critical phenomena. 
Let us assume a linear sweep of the field in a period T between Hq and hi: 

h{t) ^ ho - Vht , Vh = rp ■ 

We define the distance from the critical point, e(t) ~ h{t) — he. From the dispersion 
relation, Eq. ([3]), it follows that in the thermodynamic limit the gap in the excitation 
spectrum is 

^^\h~h,\, 

while in a finite chain with L spins at h = he the gap takes the finite minimum value 

A (or Al) sets the energy scale and the relaxation time is 

7- - Zil - '^0 



A m\ ' 

where tq is an appropriate coefficient of proportionality. The divergence of t is the 
hallmark of the critical slowing down in the usual phase transition. One can also 

define a healing length: ^ — ct, where c = lim 



= 7 is the spin wave velocity, 

h=l 



fc^TT dk 

and a density of excitations v — 1/^. 

According to ZDZ, when approaching the critical point, the dynamics of the 
system ideally should change from adiabatic to frozen behavior at time t. Since e{t) 
changes on a time scale e(t)/e{t) the crossover between an adiabatic and a sudden 
region should occur when|j 

T{i) = ^ - i ^ = To/vh . 

e{t) 

When crossing the critical point at velocity vi^ the healing length at i is ^ = ci and the 
excitations density after the sweep would be cx ,Jvh. Studying the time dependence 
of different quantities in the system we will directly probe this scenario. 



4. Effect of mesoscopic fluctuations on the ZDZ scaling law 

After a sweep from h{Q) to h{t) the density of excitations at time t in the system is 
given by the expectation value 

^=l{GS\U\t)Y,vWu{t)\GS) 

k 



X Here the time set-off is chosen such that the critical point is reached at t = 0. 
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Figure 1. Single mode excitation density vs. (pk for 7 = 1 at different quench 
rates (from bottom to top: = 0.01,0.05,0.50,1.00) and different final fields 
h( = (solid lines) and h[ = 0.9 (dashed). Dotted lines represent the predictions 
from the Landau-Zener formula. Vertical lines point out the values of {(j>k} for 
L = 129. 



k>0 

where 77^ (77) is the single mode creation (annihilation) operator as defined by Eq. ([5]) 
and U{t)\GS) = UfklHPkit) + ak{t)flf\] \Q)k,~k % At 7 = 1 and h{T) = we 
have ri^rjk — ^ '}2n{^ — cr^cf+i), so that v coincides with the kink density in the 
system at the end of the sweep |T7]. 

Fig[l] shows the excitation spectrum for different velocities. For a finite chain 
of L = 2N + 1 sites the lowest energy mode is (t>N-i = 7r(A^ — 1)/N, this means 
that, for instance, at Vh — 0.01 the sweep can be considered adiabatic up to L < 65. 
Looking at Fig[T] one can compare the exact result and the Landau-Zener formula 
Pk = exp[--27r7^ '^'"^|'^'°^ ] that reproduces surprisingly well the exact solution for final 
field h[ = even for rather high Vh- This confirms that the physics of the process is 
ruled by the Landau-Zener tunneling [ITj. One can also observe that for sufficiently 
low sweeping rate the excitation density just after the transition point {h = 0.9) is the 
same as found at /i = 0. 

An adiabatic and a sudden region can be distinguished: when the system enters 
the sudden region the number of excitations rapidly increases up to some maximum 
value. Afterwards, an oscillatory behavior of the excitation density is observed. Thus, 
a sensitive (though not unique) possibility to identify the sudden region is the time 
derivative of the excitation density ^ — —Vhj£- Using this method the sudden 
region can be depicted as the grey shaded regions in Figl^] At fixed Vh, the number 
of excitations clearly depends on the density of states (DOS) between the ground 
state energy and the "sweeping energy" £.„ :— Vh- In particular only few modes 
contribute to the total excitation density and their number is related to the minimum 
gap Al = 2tt^/L and to the DOS (upper left panel of Fig. [21). When the sweeping 
velocity increases the sudden region broadens (Fig.O right panel). The Kibble-Zurek 
scenario predict the width of this region to be proportional to the squareroot of the 

§ It is worth noticing that the fc = mode can only be singly occupied or empty. Its occupation 
corresponded to the odd sector, and therefore the A; = mode must be excluded. 
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Figure 2. Left upper panel: Density of excitations vs. = ho^^ht for different 
chain lengths (from bottom to top: L = 51, 101, 201, 401) at = 0.01 and 7 = 1. 
Dotted lines represent the prediction from the Landau-Zener formula and a finite 
number of modes. Left bottom panel: Time derivative of the excitation density. 
The shaded region points out the "sudden interval" in the sweep. Right panel: 
^ vs. h{t) for 7 = 1 and increasing sweeping velocities vj, (from bottom to top). 
The shaded areas point out the sudden regions that broadens with increasing Vf^. 




Figure 3. Left panel: for different quench rates Vf^, for 7 = 1 and L = 101. 
Full lines represent the total , and the dashed lines represent the probabilities 
to excite the first mode only. Right panel: scaling of Ah vs d/j as TV increases on 
a log- log plot. 



sweeping velocity: Ah cx v/it = y/TQVh- Already a much less detailed quantity reveals 
similar insight and admits comparison with reference jBj ; namely the total excitation 
probability 

p^ = l-\{GS{t)\U{t)\GSm\^ . (9) 

Like ly, it exhibits a maximum just after the crossing of critical point (see Fig. [3]). 
We define Ah = h^ — /imax, where /imax is the value of the field where p^ reaches 
its maximum. The value of Ah evaluated in this way depends on L due to the L- 
dependence of the gap. It is only the extrapolation of Ah to i ^ cx) that clearly 
displays the Kibble-Zurek square root behavior ^r^Vh as can be seen from the log-log 
plot in Fig. O 

In order to verify the ZDZ scaling law of the excitation density, we plot 
V as function of the sweep velocity Vh in Figdl The data extrapolated to the 
thermodynamic limit agrees well with the expected scaling law v cx ^Jvh independently 
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Figure 4. Left panel: Excitation density vs. ii;, at 7 = 1. The full circles 
represent the extrapolated A'^ — > 00 values and the dashed line is the best fit of 
those points. Right panel: the same as in the left panel, but with 7 = 0.1. 
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Figure 5. Extrapolation of the density of excitations to the thermodynamic limit 
for 7 = 1. Symbols: v{vh,L)\ dashed lines: regression lines. 

of 7. In particular for the Ising model we find v = 0.1I(4)y^i^, while previous 
analytical studies based of a LZ |17| and perturbative 13J approximated approaches 
overestimated it. With decreasing 7, the minimum gap decreases linearly in 7 while the 
density of state increases as 1 The combined effect results in an overall enhancement 
of the number of excitations in the system and in the thermodynamic limit the density 
of excitations increases as ~ I/7. This is well verified by comparison of the curves in 
FigHlfor 7=1 (left panel) and 7 = 0.1 (right panel). 

One can wonder whether the ZDZ scaling v cx ^Jvh is a pure many-body effect or 
determined already by the first excited state alone. In order to answer this question 
it is important first to observe in FiglS] the extrapolation of v for 1/L 0. The 
thermodynamic limit is reached with a perfect linear behavior in 1/L as expected by 
the constant value ^ of the DOS at low energji]]. However, this scaling holds only 
for large enough values of v^L. To be precise, a sufficiently large number of excited 
states must be taken into account in order to get the correct thermodynamic limit 
of v. Even though the first excited state certainly gives the main contribution to 
the overall excitation density at low sweep velocities, Figl6] shows that all the levels 

II The DOS in the low energy limit at h = 1 reads S^^q ~ ^ + ^^^^^ ' 
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Figure 6. Single mode contribution to tlie excitation density for L = 401 and 
7 = 1. Bold line: total excitation density; thin lines: contribution of the first 6 
modes from the lowest (in energy) to the highest one (from top to bottom). 

between the lowest excited state up to the sweeping energy = Vh are necessary for 
the ZDZ scahng to show up. At fixed vt and in the limit L oo the number of levels 
in this interval is ~ ^(^^/t ~ l-K^jV). So, we conclude the effect to be essentially due 
to many-body physics rather than to a simple level crossing. 

The last important issue of this section is to examine how the mesoscopic (finite) 
dimension of the spin chain affects the ZDZ scaling behavior and eventually the 
adiabatic passage through a quantum phase transition. In fact, two competing effects 
arise when chains of finite length L are considered: 

• on the one hand, finite L means finite gap at /i = /ic and the fidelity of the 
adiabatic passage is bounded by the ZDZ scaling law; 

• on the other hand, finite L also implies large mesoscopic fluctuations of the 
excitation number produced during the sweep, i.e. the average quantities can 
lose their importance in favor of higher moments. 

To address this question, we evaluated the mean square fluctuations of the average 
density of excitations v 

51 ^ j^{GS\UHt)J2vivkvUk' Uit)\GS) - 

kk' 

The relative fluctuations are plotted in Fig. [7] for the quantum Ising model, i.e. 
7=1, and a moderate sweeping velocity. The striking feature that emerges from this 
analysis is the high sensitivity of the relative fluctuation to the size of the system. 
Even for rather large lattice size {L = 401) the relative weight of 5^ is as large as 18% 
and rapidly increase up to 80% for L — 101. When the rate of the passage is given, 
small lattice size implies few excitations on average, but also very large fluctuations 
in their number and thus a lack of information about the final state. 

Let us remark that the ZDZ scaling emerges correctly only in the thermodynamic 
limit where the relative fluctuations are strongly suppressed. Thus, even if the 
dynamics of the model ([T]) is governed by the Landau-Zener tunneling of non- 
interacting particles, the ZDZ scaling is the proper signature of the quantum critical 
point crossed during the sweep rather than the result of single particle tunneling 
processes. 
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Figure 7. Relative fluctuations of the excitation density versus L in the Ising 
model 7=1 and sweeping velocity = 0.05. The largest lattice size is 
L = 102401. 



5. Robustness of the adiabatic passage in the presence of noise 

A real experimental setup is subject to noise, another source that may corrupt 
the adiabatic passage fidelity. Here, we examine the case of a noisy driving field, 
h^{t) — h{t) + ^(t)/2, where ^(t) is a stochastic field. The simplest choice of a 
white spectrum for i{t) (i.e. (^(t)) = and (C(i)C(O)) = T5{t)) wih be enough for 
understanding the main features due to noise and decoherence. The effect of the noise 
is that a diagonal stochastic term enters the Hamiltonian ([2]), leading to 



Hi ^ - cos 0fc IL2 



ak ^ 2 



bkCT^ . (11) 



In order to study the dynamics of such a system one has to deal with the density 
matrix formalism and to solve the Lindblad equation |25j 

Pk{t) = -t[Hr\Pk{t)]-naMa\pkm , (12) 

that can be rewritten as 

Pk = ^kPk , (13) 

where pk is the Bloch vector defined by = i(Il + jOfc -(t). The Lindbladian here takes 
the form 

/ -r -2ak 2bk \ 

Ck^ { 2ak -T . (14) 

V ~26fe y 

The stochastic noise induces a nonhermitian perturbation to the original 
Hamiltonian (jlj. The presence of such kind of perturbations can radically 
modify the physics of the unperturbed (hermitian) system, in particular when the 
latter has degeneracy or quasi-degeneracy points, as is the case here: degeneracy 
or quasi-degeneracy points typically evolve into branch-points of nonhermitian 
degeneracies [18] . Indeed, two eigenvalues (and the two eigenvectors) of the fermionic 
system (jlip merge in a finite interval of h, as compared with a single point without 
noise. This is seen in Fig. [S] where the imaginary part of eigenvalues of the Lindbladian 
Ck versus the field h are plotted. 

Analysis of the noisy case admits to answer two fundamental questions: 

• Is it still possible to discriminate the adiabatic from a sudden regime in the 
presence of noise? 

• How robust is the KZ mechanism with respect to noise? 
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Figure 8. Imaginary part of Ck eigenvalues vs. h for 7=1, fc=0.997r, and r=0.15. 
The dashed lines represent the noiseless spectrum for the same values of the 
parameters. 




Figure 9. Kink density vs. Vf^ for 7=1, L=401, and different values of the dissi- 
pative coupling: From bottom to top, T = 0, 10-*^, 10"^, 10""', 10"^, 10"^, 10'^ . 

Solving (numerically) Eq. (jl3|) we evaluate the kink density (excitation density) of 
the system using Eq. ([5]). It is shown in Fig.[n]for the Ising model (7 = 1) after a sweep 
of the field from Hq = 10 to hf = 0. The excitation density of the system increases 
considerably with the noise strength F up to the saturation value 1/2 (the maximum 
number of kinks is L/2). Besides this, there are two new and striking features: 

(i) the considerable noise instability of the adiabatic approximation even if T is orders 
of magnitude smaller than the excitation gap; 

(ii) the non monotonic behavior of the excitation density with the quench rate Vh- 

In order to understand this somehow counterintuitive behavior, all relevant time scales 
for the adiabatic passage need to be taken into account: the relaxation time t = 1/A, 
which is the only relevant scale in the noiseless case, and the decoherence time [in 
the white noise case — l/(2r)]. Both time scales have to be compared with the 
sweeping time T = (Hq — hf)/vh- In order to faithfully drag a quantum state by the 
drive, the sweeping rate has to be much slower than the gap of the system, Vh ^ A, 
but, on the other hand, it has to be much faster than the dephasing rate, Vh ^ F. 
Then, the condition for the success of the adiabatic passage is that F ^ ^ A. 

This is nicely seen from the excitation probability as a function of F (Fig. [TOl left 
panel) and explains the observed non-monotonicity. The same feature is resolved in 
terms of modes (pk in the right panel of the same figure. When F is of the order of 
Vh or even larger, most of the modes lose their quantum coherence during the sweep 
irrespective of the velocity. Eventually, the excitation density of each mode saturates 
to the value 1/2. The relevance of the appearing of the decoherence time compromises 
the distinction between the adiabatic and the sudden region. The noise progressively 
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degrades the state of the system into a mixed one. This shows up in Fig[TT] looking 
at the behavior of the purity defined as Trp^, that we find to depend solely on the 
velocity-to-noise ratio Vh/T. The rate of degradation of the state during the sweep is 
constant only far from the critical point and this is particularly evident in the case 
of low velocity-to-noise ratio (full line) where the state mixes up almost completely 
before reaching h{tc) — he- For large enough values of w/i/F, approaching the critical 
point the mixing rate of the state increase as shown by the flex-point in the purity 
close to he- This puts clear emphasis on the various facets of noise: besides the mere 
relaxation-excitation process the main effect consists in decoherence, i.e. mixing of 
the density matrix. The latter effect has been discarded in a previous work on this 
subject [19]. We conclude that the (quasi-) degeneracy point not only gives a limit 
to the rate the parameters of the Hamiltonian may change, but more importantly, it 
amplifies the effect of decoherence. 

6. Conclusions 

We have analyzed the dynamical crossing of a quantum phase transition analytically 
for an ideal clean system and numerically in the presence of white noise. Our focus 
was to determine the robustness of the adiabatic passage trough a quasi degeneracy 
point of a many-body system when one considers the effects of finite sweeping period, 
mesoscopic fluctuations and stochastic forces. 

We found a surprisingly high predictive power of the Landau-Zener formula for 
the final excitation density, when several modes are taken into account. Its validity 
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exceeds significantly the validity range attributed to it in Ref. ^17\ when the number 
of Landau-Zener modes to be taken into consideration are determined by the density 
of states and the characteristic energy scale of the sweep. This number diverges in the 
thermodynamic limit and it is only in this limit that the ZDZ-scaling [3] manifests. In 
this sense is this scaling essentially due to collective many-particle physics rather 
than to the Landau-Zener tunneling of few two-level systems. In addition huge 
fluctuations around the average excitation density dominate the state during the 
sweep, in particular for not too large systems. This tells that the physics of the 
transition is not entirely captured by the excitation average and its square root ZDZ- 
scaling. The implication of this fact deserve further analysis in order to understand its 
implications for AQC. In fact, large fluctuations on the average number of excitations 
implies a loss of information on the final state of an AQC process. 

The presence of noise corrupts the success of a possible adiabatic passage for 
low sweeping velocities due to the opening of a non-hermitian degeneracy leading to a 
"critical region" rather than just a singular point. The same phenomenon appears also 
in presence of just an avoided level crossing or a quasi-degeneracy of the spectrum. 
This leads to an all-over increase of the average excitation density, but it is shown that 
the most devastating effect of the noise on the transfer fidelity across a non-hermitian 
degeneracy region is the decay of the state purity that is emphasized by the passage 
through the degeneracy. Therefore, one has to consider a second time scale besides 
1/A: the dephasing time. These two time scales are the upper and lower bound for 
the sweeping rate. Thus there exists an optimal velocity for each noise level that 
minimizes the fidelity loss. 

Finally let us comment on different kind of noise couplings that could be of interest 
for future studies. In fact, one can wonder how the present conclusions may change 
if the stochastic variables couple with S'^(i'), breaking the symmetry of the system. 
Even though the noise coupling we considered preserves the partity symmetry of the 
Hamiltonian (but not the time reversal one), it causes both relaxation (for h ^ he) 
and decoherence. In this sense, it provides a complete phenomenology. Moreover, the 
main feature we find, namely the enhancement of the mixing up rate of the state, 
should be independent on the kind of noise coupling. With this respect, it would 
be interesting to study colored noise that can strongly modify the dephasing time; it 
should also have an effect on the white-noise induced increase of the excitation density 
with time during all the sweep. 
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Appendix A. Exact diagonalization of the XY model in transverse field 

In order to diagonalize the one-dimensional XY model in transverse magnetic field ^ 
we follow Refs. [71 [M] [15] and perform the substitution 



a] :=S^ + ^Sy 



(A.l) 




(A.2) 
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The resulting (hard-core bosonic) Hamiltonian becomes 



^h{t) n, + \Lh{t) .{A.?,) 



L=2N+1 L=2N+l 

^ (a]aj+i + h-c.) + + h.c.) 

j=i J j=i 

Here we restricted ourselves to an odd number of sites L = 2N + 1 with integer 
N. Though this is not crucial for the exact diagonalization of the Hamiltonian, this 
becomes equivalent to the restriction to an odd total number of particles - and hence 
no boundary phase - when starting from a fully polarized state at very large magnetic 
field. 

The Jordan- Wigner transformation 

/ j-i \ 



a'- exp ITT 



1=1 



further transmutes the Hamiltonian into fermionic form 

L=2N+1 



H = — 



1 



L=2N+l 



Subsequent Fourier transformation 



(A.4) 
(A.5) 

-Lh{t) (A.6) 
(A.7) 



completely decouples the Hamiltonian, so that H = ®^^.^Hk, where Hk acts in the 4- 
dimensional Hilbert space {|0)^, , \k, —k) ; \k) , |— fc)}. Due to the parity symmetry, 
it further decouples into Hk = Hf'^ © Where Hf'^ = - cosffiA, while in the 

even-occupation Hilbert space and basis {|0) , — fc)} one obtains 

h{t) -Z7sin0fe 
i^ sin (f)k — (2 cos 0^ + h{t)) 

cos -t- h(t) —i"f sin (p^ 
i7sin(/)fc — {cos (j)k + h{t)) 

Qk -ibk 
ibk -Ok 



COS (t>k I2 



= - COS (f)k I2 

Expressed in Pauli matrices this is 



(A.; 



jjeyen 



COS (t)k^2 + OfcCT^ -I- bkCr'^ 



The eigenvalues and the ground state of the nontrivial second part are given by Eqs.Q 
and ([3]), respectively. 



Appendix B. Exact time evolution of the XY model in transverse field 

The Heisenberg equation ([7]) for the time evolution operator has been solved in Ref. [7] 
to obtain the exact time evolution of correlation functions of the model. Using the 
notation 

/ Tj{k) Tj{k) \ 
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we obtain the following set of differential equations 
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(B.l) 
(B.2) 



^From these equations and the initial conditions one concludes C/22 



U. 



(k) 



21 



11^2^* and the missing amplitudes are extracted as 



U- 



(fc) 
12 



bk 
1 



(fe) 



22 iO'kU 22 



Ull>' and 

(B.3) 
(B.4) 



The differential equation (jB.ip is exactly solvable if Ok is polynomial in t [23] and 
even for particular cases of transcendental dependence of t [7] . For the simplest case 
h(t) = ho — Vht, with Vh — {ho — hi)/T, the solution is 



^11 — e 



4wh 



/2^ 



{h{t) + cos (t>k)] 



+C2F[ir 



2Cos2(?!)fe 



■^—^,^,i2vh{h{t)+cos^k)f] 



2COS20fe - 1 

C^iH «7 

4uh 



1 + i 

1, +COS0fc)] 



/2^ 



+C2F[zf H^i^^ - i i z2^;,(/.(t) +cos0,)2] 

In our notation F[x^y, z\ is the confluent hypergeometric function and has the series 
expansion 

{x)k 



F[x,y,z\ = ^ 



fe=0 



{y)k kl 



with {a)k '■— a{a + 1) • • • (a + fc) the Pochhammer symbol. Its pole structure is that of 
the gamma function r(y) such that F[x^ y, z] — T[y]F[x, y, z] with the entire function 
F[x,y,z]. The generalized Hermite polynomials H^[z] =: H[iy,z] are expressed in 
terms of the confluent hypergeometric function 



H[v,z] = 2"^/^ 



Fj 2 ' 2 ' ^ ] 



P[ 2 ' 2 ' ^ ] 



At integer values of ly, this function coincides with the Hermite polynomials. 
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